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The Rosensweig instability has a special character among the frequently discussed insta- 
bilities. One distinct property is the necessary presence of a deformable surface, and another 
very important fact is, that the driving force acts purely via the surface and shows no bulk 
effect. These properties make it rather difflcult to give a correct weakly nonlinear analysis. 
In this paper we give a detailed derivation of the appropriate amplitude equation based on 
the hydrodynamic equations emphasizing the conceptually new procedures necessary to deal 
with the distinct properties mentioned above. First the deformable surface requires a fully 
dynamic treatment of the instability and the observed stationary case can be interpreted as 
the limiting case of a frozen-in characteristic mode. Second, the fact that the driving force is 
manifest in the boundary conditions, only, requires a considerable change in the formalism 
of weakly nonlinear bifurcation theory. To obtain the amplitude equations a combination of 
solubility conditions and (normal stress) boundary conditions has to be invoked in all orders 
of the expansions. 

Subject Index: 034, 054, 056 

§1. Introduction 

Since its discovery in 1967, the normal field or Rosensweig instabilitjl3 attracted 
the attention of experimentalists and theorists, alike. The phenomenon describes the 
transition of an initially flat ferrofluid surface to hexagonally ordered surface spikes 
as soon as an applied magnetic field exceeds a certain critical value. Ferrofluids are 
suspensions of magnetic nanoparticles in a suitable carrier liquid. They are coated 
by polymers or charged in order to prevent coagulation and show various distinct 
material properties o One of the most prominent examples of these properties is 
the superparamagnetic behavior in external magnetic fields, which accounts for the 
large magnetic susceptibility and the high saturation magnetization in rather low 
magnetic fields. If one starts the cross- linking process in a mixture of a ferrofluid and 
a polymer solution with cross-linking agents, a superparamagnetic elastic medium, 
called ferrogel, is obtained!^ As in usual ferrofluids, the initially flat surface of 
ferrogels becomes unstable beyond a critical magnetic fieldS 

With its discovery, a first theoretical description together with a linear stability 
analysis was given0 At the free surface the stabilizing forces of gravity and sur- 
face tension compete with the destabilizing magnetic force. Although the applied 
magnetic field is homogeneous and therefore no net-force is acting on the medium, 
fluctuations of the surface lead to focusing effects rendering the local fleld at the sur- 
face inhomogeneous. With this model, the prediction of the critical magnetic field 



typeset using VTpT^.ds (Ver.0.9) 
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and the characteristic wavelength, which turned out to be the capiUary wavelength, 
was possible. The linear growth behavior was discussed later onoTlP For magnetic 
gels the elastic force contributes as a stabilizing effect leading to a shift (as a func- 
tion of the elastic shear modulus) to higher magnetic fields whereas the characteristic 
wavelength remains unchanged,^' First experiments to confirm the threshold shift 
are performed using thermoreversible magnetic gel^ as the magnetic medium. 

A linear stability analysis provides us just with the threshold and the most 
unstable mode. No prediction of the arising pattern can be made, nor does it give 
the dynamic behavior beyond the threshold. A nonlinear analysis of the Rosensweig 
instability, however, turned out to be very complicated mainly due to the fact that 
the instability necessarily involves a deformable surface. In 1977 a very first approach 
to the nonlinear regime was given by A. GaihtisI?! Since the pattern formed is 
static, Gailitis discussed the surface energy density, consisting of the gravitational 
energy, the energy contribution due the surface tension and the magnetic energy, as 
a function of the deflection of the surface from its flat ground state. Prescribing the 
regular surface patterns of stripes, squares and hexagons, he found upon minimizing 
the energy density, that at the linear threshold hexagons are the stable configuration 
that in turn transform into squares upon further increase of the magnetic field. Both 
transitions are accompanied by hysteretic regions. The stripe configuration instead 
is always unstable with respect to one of the other two patterns. A major drawback 
of this method is, that it is valid only in the asymptotic limit of vanishing magnetic 
susceptibility. Friedrichs and Engell2l extended Gailitis' method to systems with a 
finite depth and additionally gave an estimate of the maximal magnetic susceptibility 
up to which the method gives reasonable results. The energy method was extended 
to the Rosensweig instability in isotropic magnetic gels considering additionally the 
elastic surface energy densityEl} Another unsatisfying aspect of this method rests in 
the fact, that it cannot describe growth rates, since it ignores dissipative processes. 

Based on the linear result of a static instability, an ex pan sion of the basic static 
equations governing the ferrofiuid behavior was discussed EH' The analysis is, in 
contrast to the energy method, valid for any given magnetic susceptibility, however, 
for typical values no stable pattern could be found at the linear onset. Another ap- 
proaclllil' considers a static regime, where only the normal stress boundary condition 
is considered for the nonlinear expansion and where a horizontal field component of 
the magnetic field was assumed to be strong enough to suppress two dimensional pat- 
ternsii^ The dynamics of the system has first been taken into account by Kubstrup 
et alJi^ who used a Swift-Hohenberg model to describe fronts between hexagons 
and squares. This approach, however, lacks the connection of the parameters intro- 
duced in the Swift-Hohenberg equation to the material properties of the medium. 
In addition, it is unclear whether the terms involving time derivatives in this model 
are the appropriate ones. 

What one would like to have is a systematic nonlinear expansion of the basic 
hydrodynamic equations in analogy to Ref. 17). To adapt this method in the case 
of the Rosensweig instability, the adjoint linear eigenvectors in the presence of a 
deformable surface are needed to satisfy Fredholm's theorem. To circumvent Fred- 
holm's theorem, Malik and SinghEH'llll' restricted their discussions to potential flows 
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only. However, as can be seen in the present paper and in Refs.|4|), 11 ), the rotational 
flow contributions are needed to guarantee the free surface to be stress- free. Recently, 
the adjoint system for the Rosensweig instability in isotropic magnetic fluids and gels 
(and for the Marangoni instability) in presence of deformable surfaces was given by 
the present author3H2l as a prerequisite to access the nonlinear regime via a weakly 
nonlinear analysis. The latter, however, cannot be carried out straightforwardly, 
since the driving force acts on the surface, only. As a result, the (bulk) solvability 
conditions are not sufficient to give the amplitude equation, but have to be combined 
with the (normal stress) boundary conditions. This new procedure will be discussed 
in detail and executed explicitly in this manuscript, where we concentrate on the 
case of magnetic gels. 

Deformations of the free surface are crucial for the Rosensweig instability. The 
kinematic boundary condition relates the temporal changes of the surface deflection 
with the velocity of the bulk material at the surface. An a priori static description 
would, thus, completely miss this important boundary condition - one of the reasons, 
why previous attempts to deal systematically with the nonlinear instability regime 
have failed. For gels, in addition, the dynamic coupling between elastic deformations 
and flow also requires a fully dynamical treatment, even if the final instability is 
stationary. Only at the end the static limit can be performed. The pattern observed 
by the Rosensweig instability is thereby characterized as the limiting case of a frozen- 
in surface wave mode. 

Our discussion is organized as follows. In §§2 and 3 we introduce the basic hy- 
drodynamic equations and their general expansion into the nonlinear regime. Special 
emphasis is put on explaining the consequences of Fredholm's theorem for the present 
systems. In §4 we solve the second perturbative order. The necessary solutions in 
the third order are derived in §5, whereas in §6 we give the amplitude equation for 
the Rosensweig instability in isotropic magnetic gels followed by a short discussion 
of the amplitude equation for magnetic fluids in §7. Many of the detailed algebraic 
calculations are put into appendices. Some of the results have recently been pre- 
sented in a conference proceeding^lll without, however, laying out their derivation 
and the non-standard subtleties involved. 

§2. Basic equations and general approach 

The basic equations we are concerned with when discussing the Rosensweig in- 
stability in magnetic gels, are the hydrodynamic bulk equations derived for isotropic 
magnetic gels by Jarkova et alB together with some approximations discussed 
below. 



dtgi + d,T,j = pGi (2-1) 
{dt + Vkdk)eij - ^ {diVj + djVi) = (2-2) 

^^Vi = (2-3) 



diBi = (2-4) 
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^.i.kdjHk = (2-5) 



They account for the conservation of hnear momentum (2-1), the conservation of 



mass (2-3) and the fact that the elastic network breaks continuous translational 



symmetry (2-2). In our notation g is the momentum density, v the velocity, p the 
pressure, G represents the acceleration due to gravity and B and H are the magnetic 
induction and the magnetic field, respectively. The second rank tensor €ij denotes 
the strain field, while the material parameters /i2 and f2 stand for the shear elasticity 
and shear viscosity, respectively, and are contained in the stress tensor T. 

The underlying assumptions are as follows. Even though the magnetic field is 
considered a slowly relaxing variable in the hydrodynamic theory of Jarkova et al., we 
assume that it relaxes fast enough on the time scale considered in our discussion of the 
Rosensweig instability. This is justified by the fact, that the growth of surface spikes 
takes place on a time scale long compared to the temporal variations of the magnetic 
field. The magnetic field is then defined by the static Maxwell equations ( 2-4|2-5 ) 



and the corresponding boundary conditions at the surface. We also assume, that 
the macroscopic material parameters like the shear modulus and the shear viscosity 
are independent of the magnetization in the medium. This also implies that we will 
neglect magnetostriction in our discussions. Furthermore we assume the magnetic 
gel to be incompressible {p = const., en = 0) and be described by linear elasticity 
theory. Although a realistic quantitative treatment of polymeric gels requires the 
use of nonlinear elasticity, there is no reason to expect the elastic nonlinearities to 
change the qualitative behavior of the Rosensweig instability. 

The stress tensor of the magnetic medium is defined via the conservation equa- 



tion for the momentum density (2-1) and is given in our notation by 



Tij — QiVj + p6ij — {BiHj — -BkHk 6ij) — fJ,2{^jk^ki + ^ik^kj) — '^P-2^ij 

-V2{djVi + diVj), (2-6) 

while the vacuum stresses are solely due to the magnetic field for what reason the 
stress tensor there reduces to the known vacuum Maxwell stress tensor T^'^iHl' 

The hydrodynamic and magnetic bulk equations are supplemented by boundary 
conditions at the deformable surface defined by z = ^. Aside from the usual magnetic 
boundary conditions the tangential components of the mechanical stress between the 
magnetic medium and the vacuum above is required to vanish at z = ^, while the 
normal stress difference is balanced by gravity and surface tension. 

n X T • n = n X T™^ • n (2-7) 

n • T • n - n • T^^'^ • n = ordivn - pG^ (2-8) 

n X H = n X H™^ (2-9) 

n-B = n-B™^ (2-10) 

where H™'^ and B™^ denote the magnetic field and the magnetic flux density in the 
vacuum, respectively, gt is the surface tension, and where we introduced the surface 
normal n = d{z — '^)/ | d{z — Additionally, due to the deformable surface, we 
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Fig. 1. Small periodic perturbations ^{x,y,t) of the initially flat surface z — between the fer- 
rogel/ferrofluid of susceptibility x in the lower half space and vacuum in the upper half space. 
The magnetic field is oriented parallel to the 2— axis, whereas gravity is pointing in the opposite 
direction. 

have to consider the kinematic boundary condition modehng the dynamics of the 
free surface at z = ^ 

dtC = V,. (2-11) 

Performing a weakly nonlinear analysis of the stationary state evolving slightly 
beyond the linear threshold Mc, we have to expand the macroscopic variables in 
terms of e, the normalized difference of the actual applied magnetic field to the 
critical one 

{p, B, H, M} = {po. Be, He, Me} + e{p^'\ B^^\ U^^\ M^^^ + ... (2-12) 
{v,e,„e} = + e{vW,eg\e(i)} + ... (2-13) 

The magnetic field, however, is an externally given parameter acting as the control 
parameter, the series expansion of H can therefore be reinterpreted as the definition 
of e. The linear threshold is giveilll' by 

m2 = i±i^ (ly^cJTpG + 2//2) , (2-14) 

which also determines the critical fields Be, He via the linear magnetic constitutive 
equation B = H + M = /xH employed here. 

In our linear discussionjll the surface deflection ^(x, y, t) was modeled using 
plane waves ^{x,y,t) = ^g*^*-*k-r ^ nonlinear discussion, this ansatz has to be 
expanded. The linear description can just provide the characteristic mode becoming 
unstable at the threshold. The most general ansatz as a starting point for a nonlinear 
discussion is to assume of these characteristic modes with diff'erent orientations. 
Each of these modes i consists of a right and left traveling contribution (subscripts 
R and L, respectively) 

N N 

i i 
N 

= J2^iiRe''^'^~'^'''' + iiLe-'^^^-'^^-"" + |*jje-''^'*+^'^'-'" + ^"*^e^^'*+^^»-'"), (2-15) 
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Fig. 2. The sketch shows the relative orientation of the wave vectors under consideration in the 
amplitude equations ( 6-17|6-19 l. It allows to discuss the stability of hexagons and squares and 
their interaction. 



where the asterisk denotes the complex conjugate and kj characterizes the direction 
of the i-th mode. The wave number k = |k| = |kj| is the same for aU modes. At the 
end it will turn out that basically three patterns are important, hexagons, squares 
and rolls (or stripes). They are described by six critical wave vectors, for which 
we choose the geometry of Fig. [2j This geometry allows us to discuss hexagons 
(6 =6 = 6/0 and ^4 = ^5 = ^6 = 0), squares (6 = 6 / and ^ = for 
i G {2, 3, 4, 6}) and rolls (6 / and 6 = for i G {2, 3, 4, 5, 6}). As discussed already 
in the derivation of the adjoint systemp^ we have to treat the system dynamically 
and perform the limit towards a static system in the very end only. The eigenvectors 
in linear order are known to be modulated by ^(^113 and therefore separate into left 
and right traveling contributions together with the complex conjugates, similarly as 
in Eq. (2T5). The corresponding coefficients depend on the vertical direction ~e'^^ 
and ~e'^^ and are given in Ref. 11). 

While performing a weakly nonlinear analysis, we have to specify the scales 
in space and time. In a first approach we will assume a surface pattern that arises 
homogeneously in space, which allows us not to rescale the spatial degrees of freedom. 
Time, however, will be rescaled in the following manner 



t 



(1) 



et and t 



(2) 



(2-16) 



which will lead to the substitution for the time derivative 

a^^a(o) + ea« + 62ar^ + ... (2-17) 

We can interpret the scaling in time in the sense, that the dynamics of the amplitudes 
itself takes place on the slower time scales, 6,{i?,L} ~^ 6,{_r,l}(^^^M^^^ • • • )• 
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§3. Fredholm's theorem and the adjoint system 

With the scaling of time and the expansion of the macroscopic variables in terms 
of e that we have introduced so far, the whole system of differential equations can be 
expanded in terms of e. Let Cq be the linear operator and | ^) =\^^^'^) + e \ 'F^^^) + . . . 
the macroscopic state vector. The different orders in e are then given successively 
by 

£0 = (3-1) 

£0 lO^^^)) = |AA(^'W,-P'«))+ |r(5f^i^'(^))) (3-2) 



The first equation ( |3-l[ ) represents the linearized set of equation used in the discussion 
regarding the linear stability HJ Furthermore Eq. (3-1 ) defines the kernel of the linear 



operator £0, given by the linear eigenvectors | ^^^^). In the second perturbative order 



the set of equations (3-2) becomes inhomogeneous due to the nonlinear nature of 
the basic set of equations (represented by •)) and due to the rescaling of time 
(represented by ?'(•)). In the case that these inhomogeneities reproduce elements of 



the kernel of the linear operator Cq, equation (3-2) cannot be solved. The requirement 
that the inhomogeneities have to be orthogonal to the subspace spanned by the linear 
eigenvectors | ^) provides us with an additional solvability condition. It is named 
after Fredholm and reads in the second order 

A^(<f(i),<P'(i))) + {^\ r(af^^'(^))) = 0, (3-3) 

where (o | b) denotes the suitable scalar product 

L L C T 
(a|6)=^lim ^ jdx jdy j dz jdt a* b (34) 

-L -L -00 

taken over the whole range of the deformed sample. Application of this scalar product 
requires to explicitly expand all boundary values in terms of the surface deflection 

The derivation of the required adjoint eigenvectors \ for instabilities with a 



deformable surface is given in Ref. 20). Here we recall the results needed for the 



upcoming calculations. The components of the adjoint velocity field are given by 



(3-7) 



q^ + k'^ J q^ — k"^ 

where the adjoint frequency a) is given by —oj and the adjoint inverse transverse 
decay length q is defined by cf = k'^ — puP' /{^2 — ^^1^2) in the same way as the non- 
adjoint one, q^ = k'^ — pto'^ /{fj,2 + iu}i^2)- Here, ki^x and ki^y are the x and y component 
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of the wave vector of mode i. The corresponding components of the adjoint strain 
field turn out to be 

= (3.8) 



... = -4M.-t^(."-^.-Jp3p«. (3.11) 

—2 7 2 

ey, = -iiti2h,y(^e''' - e^"^) ^^^^ ■ (3-13) 

§4. The second order 

The fact that within our assumptions the magnetic bulk equations completely 
decouple from the hydrodynamic bulk equations, has two important consequences. 
On the one hand this allows us to discuss and solve these two systems subsequently, 
i.e. we first solve the magnetic part in a given perturbative order for a given surface 
deflection ^, and feed back this solution into the respective order of the hydrody- 
namic system. On the other hand, however, we have to face the problem, that 
the control parameter (the magnetization or the magnetic field in our case) does 
not occur in the hydrodynamic bulk equations and that the bulk equations for the 
magnetic system are homogeneous in all perturbative orders, which makes it impos- 
sible to obtain the control parameter in the next order by Fredholm's alternative, 
only. The coupling between these two systems is, however, mediated by the surface, 
and more precisely by the normal stress boundary condition. Satisfying the normal 
stress boundary condition provides us with an additional condition supplementing 
Fredholm's theorem. 

4.1. Magnetic contributions 

We start by solving the magnetic system of bulk equations together with the 
corresponding boundary conditions. The external magnetic field is the control pa- 



rameter defining the expansion (2-12). Due to the deformations of the surface the 
actual magnetic field will be subject to perturbations and we separate the total mag- 
netic field into the applied field H and the distortion field h. These perturbations 
still obey the linear electromagnetic equations (b = and d h = = dxh), which 
allows for the introduction of a magnetic scalar potential h = —d^. The perturba- 



tion potentials 'P (and in the vacuum) are also expanded according to Eq. (2T2) 
and fulfill the Laplace equation. In second order we obtain 

Z\<?(2) = and Z\^(2)vac ^ Q (4.^) 
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in the medium and in vacuum, respectively. In the second order of the e— expansion 
the magnetic boundary conditions for the tangential component of the total magnetic 
field H + h can be simplified to (cf. App. |A]) 

a^<j(2)vac _ Q^^(2) ^ J^M, 9,^(2) ^ (4-2) 

1 H~ 
1 H~ 

while the boundary condition for the normal component of the flux density B reads 

where we introduced abbreviations that depend on the angle Oij between the i-th 
and the j-th main characteristic mode 

kiij = + 2 cos Oij (4-5) 

k2ij = kc\/2 — 2 cos Oij . (4-6) 

A convenient ansatz for the magnetic scalar potentials to solve this system of 
equations consists of two contributions. The first contribution ^(^,1) proportional 
to the linear deflection ^^^^ to account for the contributions proportional to M^^^ in 
the boun dary conditions (4-2) and (4-3). This automatically satisfies the Laplace 
equation (4-1) for ^(2,1) j4| fj^j^g second contribution accounts for the higher 

harmonic couplings of the linear characteristic modes proportional to ^^'^\ which are 
modeled by the product of two characteristic modes 

i^^^ =kcY.{i^^,+^^^*+C.C.). (4-7) 



The characteristic wave vector in Eq. (4-7) is just added to give the same unit 

as e^. 

The Laplace equation (4T) for (?(^'^) is satisfied by the ansatz 



^(2,2) = fc, ^(^g-'^^Cieie'^"^- + ^^^^^""iii^e^^^^' + c.c.) (4-8) 



and by a corresponding one for the magnetic potential in vacuum. 

The boundary conditions for the different Fourier modes decouple and can be 
satisfied separately. We obtain for the contributions proportional to ^'•^^ 

^(2,1) ^ _i^^(l)gfcc^ (4.9) 
1 + /U 

(^{2,l)vac _ 

which are of the same structure as in the linear case. The presence of M(i) guarantees 
to be of second order. 
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The contributions due to the higher harmonics of the characteristic modes read 



(2,2)_R 



^(2,2)iivac 
ij 



(1 + /^)^ 
,2 



Mr 



lij 



(4-11) 
(4-12) 



while ^jj'^''^ and ^j^'^^^™'^ are obtained replacing kuj by k2ij in Eqs. (4-11 ) and Eqs. 
(4-12), respectively. 



4.2. Hydrodynamic contributions 



According to the general expression ( 3-2 ) the set of hydrodynamic bulk equations 



is given in the second perturbative order by 



pd^'^^vP + _ 2f,2dje\f - U2{djdivf + 8^8. 



,(2) 



-pdi^v^^-dApv^Kf^-2,2^^) (4.13) 



d?^4?-lid.vf + 8,vr) = -8^^e^-v^^8,4^ 



a (2) 



0. 



(4-14) 
(4-15) 



The structure of these equations suggests two kind of solutions similar to the mag- 
netic part. One contribution is proportional to the main charac teris tic modes S,^^^ 
and a second one proportional to the second harmonics ^^2)^ gq, (4.7). 

The corresponding boundary conditions at the surface z = are expanded in the 
same manner (for a detailed discussion cf. App. |B]). For the tangential contributions 
we obtain 



2/i26(2) + ,.2{8,vP + 8,vP) = 12(2) , 



(2) 

where the inhomogeneities are abbreviated by and are listed in App. 
( B 2 ) and ( B 1 ) . In contrast to the gravitational and the magnetic force the' 



(4-16) 
(4-17) 

Eqs. 

surface 



B 



tension always acts normal to the surface and therefore enters only the normal stress 
boundary condition 



1 + /i 



(4-18) 



(2) 

with Qzz given in Eq. (B 3). Finally, the kinematic boundary condition describing 
explicitly the deformable surface reads in second order 

9(0)^(2) + 0(1)^(1) + . 9) ^(1) = yf) + £_i^)8jp . (4-19) 



The last contribution in Eq. (4-19) is due to the fact, that in second order the surface, 
at which the boundary conditions have to be evaluated, is already deflected. 
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4.2.1. The solvability condition in second order 

The general solvability condition discussed in §A[3]is applied to the set of second 
order equations (4T3 4T5) and explicitly reads 

{Vi I -pd^- - dj{pv- V- - 2p2ej,, e^. j) 

+{e^J\-^l'\f-vl}^^ke^^)=0. (4-20) 

At this point one might be tempted to use the fact that the Rosensweig instability is 
a stationary one (in linear approximation) and substitute = as well as 

the stationary limits of the adjoint and original eigenvectors into condition (4-20). 
The solvability condition would then reduce to 

i^ij I ^4'^) = (± ^^^'^ + ^^'^^('^^ I ^^^■) = (4-21) 

corresponding to the solution lo^^^ = = a^^\ Here, we have replaced d^^^ by 
±iuj^^^ + a^^^ (for right- and left-traveling waves, respectively) implying a normal 
mode ansatz for the time dependence of the amplitudes. Of course, u;^"-* = is 
the correct solution in the stationary limit. However, in that limit the connection 
between bulk equations and boundary conditions is lost (cf. Eqs. (2-2) and (2-11)) 
and an amplitude equation cannot be derived. Therefore, one must still treat the 
system as fully dynamic at least at those places related to the kinematic boundary 
condition and to the velocity/strain relation, and satisfy Fredholm's theorem with 
the derivative dj:^^ being finite. One can, however, at non-crucial instances simplify 
the calculations by the fact that ^^(0) is small, but only at the very end one can take 
cuW = 0. 

The solvability condition (4-20) consists of two different parts. One containing 

spatial derivatives and the other the (scaled) time derivative We first discuss 
the latter part. The integration upon x and y is straightforwardly done and only 
retains contributions that are proportional to 5(kj — kj). After integration with 
respect to z we end up with the following expression. 



lOJ 



U t ^2iuit t U -2iu}t\ ^2at ( o , , ^c(fcc + g^)^ 

^^d^Re -U^^ne )e {8p2 ^(^^ + 

"^^f"' ^ "f" ^ ^ Wk^TW J 



.r,(0)l2 r (0)12^ 4feg + Gk^q + 6ktq^ + 6/cgg3 + afcgg^ 

-^^["^ ^ ^ ^ qkl{k^ + qr / • ^^'"'^ 

For the second order contributions we finally get 

{v, I 9f)(p^;f))) + (e-|af)e;])) 
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+a^'hf^2k,{iU,R + US*R + LdlR + Ldk) , (4-23) 

where the static Umit has safely been performed. 

Up to now it has been possible to do the calculations without specifying the 
actual number of modes contributing to the nonlinear pattern and the results are 
applicable for any value of and in particular for any angle between these modes. 



This is changed when the second part of Eq. (4-20)), containing the spatial deriva- 
tives, is considered. Two of these terms turn out to be irrelevant for the second 
order solvability condition and are not discussed here. The only relevant term, 
2fi2{vi I djie^jl^e^j!'^ )), generally vanishes, except when three linear modes oriented at 
it/ 3 relative to each other are interacting. This hexagonal order is enforced by the 
integration upon x and y. Integrating with respect to z yields in lowest order of uj^^^ 
and ct(0) 

-?1l6h?3L + 6r6l'?3L - C.C.) 
-3(T*^°V2^c (6-r|2_R^3/? + Iil|2L^3L + Ii_r|2k6l 
+CiR^2l£,3R + £,1l£,2rC3R + S,1l£,2lC3R 

+llL6Ke3L + ilRLLi3L + C.C.) . (4-24) 



Eqs. (4-23) and (4-24) are the two parts that enter the solvability condition Eq. (4-20), 
which we are now going to solve. The imaginary part yields the condition 

4iw(i) iikiiR-iiLim) = -3iw(°) A:c(liij6ije3ij-eiLl2Le3L+eii?6Re3L 

+^1R^2L^3R - £,1L^2rC3R - ^1lC2lC3R 

-ilL^2Ri3L + ilRi2d3L-C.C.) . (4-25) 

This condition is identically fulfilled by the ansatz 

iiL = iR = L and i*^ = i*j, = i*, (4-26) 

which is the solution one expects for the stationary case, since in that limit one 
cannot distinguish right from left traveling waves. 

Using this result for evaluating the real part, we obtain 

2a(i) ^^^^ = -S'^^^^^cia^^ + iim) , (4-27) 

i 

which obviously is solved by 

a(.^)^, = -a^o)kJ*i* and llip = l^sP = II3P (4-28) 
and all its cyclic permutations 1 — )• 2 — )• 3 — )• 1 and their complex conjugates. Equa- 



tion (4-28) tells us, that the slow variable a^^-' scales in the bulk with a^^\ indicating 
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that a^^^/a^^^ stays finite in the stationary Hmit. This behavior is mediated by the 
the kinematic boundary condition dtC = Vz (2-11). As a consequence, the velocity 



field as well as the adjoint velocity field are proportional to the time derivative as we 



realized in Ref. 11) and in Eqs. (3-5)-(3-7). This is physically reasonable, since in the 



case of the Rosensweig instability the velocity field vanishes if the surface pattern has 
fully developed and the hydrodynamic bulk equations are trivially fulfilled by v = 0, 
the same solution as for the initial undeformed ground state. This singular behavior, 
unique for the Rosensweig instability, is scaled out by the choice of a dimensionless 
time derivative dj:^^ = a^^^ /a^^^ for the bulk hydrodynamic equations. Using this 



time derivative, Eq. (4-28) can be rewritten as 



(4-29) 



Equation (4-29) gives the relation among the three amplitudes of the second 



order deflection, characteristic for hexagon patterns. For any other regular 



pattern the right hand side of Eq. (4-24) is zero implying, that there is no nonlinear 



interaction between two different modes in the second order for those patterns. 



What is missing in Eq. (4-29), which in a sense can be viewed as a primitive 



form of an amplitude equation, is a contribution proportional to the control pa- 
rameter M^-*^). This is due to the fact, that the two bulk systems of magnetic and 
hydrodynamic equation decouple completely. The control parameter enters the am- 
plitude equation via the normal stress boundary condition, the only way magnetic 
and hydrodynamic subsystems are interacting. 

4.2.2. Solutions proportional to ^^^^ 

we have to 



Before we can exploit the normal stress boundary condition in ^4.3 



determine the solution of the hydrodynamic contributions, Eqs. (4-13) - (4T7). From 



Fredholm's theorem we learned, under what conditions we can flnd a solution to the 
system of equations in the second perturbative order. As in the magnetic part, we 
distinguish solutions of the system of equations that are either proportional to ^'•^^ 
or proportional to ^^^^ . In this subsection we concentrate on the part proportional to 
Inspired by the linear discussion, we use a scalar (^(^'^) and a vector potential 
ijr(2,i) fQj, ^]^g potential and the vorticity flow, respectively. For the contributions 
proportional to the main characteristic modes the governing equations read 



(0),7v(2,l) 



,(2,1) 



(0)^2^(1)^(1) 



(4-30) 
(4-31) 
(4-32) 



with the abbreviation /i2 = M2 + 1^2 <9j^'*. On the right hand side of these equations 
the first order (linear) potentials act as inhomogeneities. They are listed in Ref. |4]). 
The appropriate boundary conditions for the fiow potentials are derived in App. 



B.2 and read for tangential stress 



(2,1) 



ji2{dl - 5^)0^(2,1) ^ ^,9^9^^(2,1) _ 2/i2a.a.^(2,i) 




0. 



(4-33) 
(4-34) 
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The physical boundary conditions have to be taken at z = ^^^^ in the second order. 
This leads to additional contributions in which have already been taken into 
account in the effective boundary conditions Eqs. (4-33) and (4-34). The latter 
therefore have to be taken at z = 0. 

The kinematic boundary condition now involves the slow timescale t^^^ and reads 

^(2,1) ^ 5(1)^(1) , (4.35) 



We start with the particular inhomogeneous solutions of Eqs. (4-31) and (4-32) 
for the vector potential ^ and the pressure p, respectively. It can be checked that 
the following fields satisfy the inhomogeneous bulk equations 



(2,1) ^ ^(2,l)inhom^(l)^^g2 



C^^he''^ and p(2.1)inhom ^ _^Q^ 



(4-36) 



with the operators defined by 

^(2,l)inhom^_^^2^^(l)^^ ^(2,l)inhom ^ M2^g(l)^^ ^ (4.3^) 



/^23 



The i nhomogeneous solutions do not y et sat isfy the boundary conditions (4-33) and 
(4-34). Substituting *inhom jj^^^ (4.33) results in an additional source of tan- 
gential stress at the boundary due to the inhomogeneous solutions, which can be 
balanced by the homogeneous ones 



^2 



(4-38) 



If we use the following ansatz for the homogeneous solutions of the flow potentials 

^(2,l)hom ^(2,1) 



^(2,l)ho 



^^2,l)hom ^ g^^(2,l)gg.^(l) 



and ^ ^{2,l)g/c..^(l) 



the amplitudes tf'^^'^) are given by 
^(2,1) ^^'^ 



^(2,1) _ A^2 + ^2 g(l) 



(4-39) 



(4-40) 



Note that q is the inverse decay length of the linear transverse modes with = 
k1 + p{df^]^ /{^i2 + 1^29^^^ and S^^-* is a short hand notation for ztiw^-^^ + (t'^^\ as 
before. 

The homogeneous solution of the pressure p(2:i)hom jg straightforwardly given by 
Eq. (f^sTl) 



P 



(2,l)hom 



(4-41) 
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and if we exploit the kinematic boundary condition (4-35), the solution of the scalar 
flow potential ip^"^'^^ can be determined as 



(2,1) 



(1) 



2k 



kc{q^-k^,) V * 'f,2{q^ + k^,) * 

With the help of the flow potentials the velocity fields are determined 



(4-42) 



,(2,1) 



2 2^^2+/f2 ,2 
^2 



kcz , n/^2,2 „2 /^2 + /i2,2/ 2 



Ai2 fJ'2 



kt{q-ki)- 



q J q^ — kl 



,(2,1) 



Mq^-k'i) 

with the abbreviation 



L(-)a,"'{.'" and if." 



fi2[q^-klY 



(4-43) 
(4-44) 



L(^) =[/i2(g'-/c2)-2/X2/c2 



+ [2fiq^ - (/i2 - /i2)(g' - + qz)] 



(4-45) 



from which the strain fields follow 



.(2,1) 



(2,1) _ fJ'2 + A2 



ki,akifiL-{z) J„^'^J■ 



(1) 



/^2 

/"2 + A2 



(4-46) 
(4-47) 



" 2/i2 U2-/t2 



(4-48) 



for {a, G {j;, y} with 



q^ — k'l 



k.e^^'-qe'^' 



^l + qz^g. 



(4-49) 



This concludes the derivation of the second order eigenfunctions that are pro- 
portional to These solutions satisfy every condition except the normal stress 
boundary condition. The latter will be used to determine the still unknown first 
order correction to the control parameter, M^^\ which finally enters the amplitude 
equation as the linear contribution. We postpone the actual derivation of these 
contributions to ^4.31 
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4.2.3. Solutions proportional to ^^^^ 

We are left with solving the system of hydro dynamic equations in the second 



perturbative order, Eqs. (4-13 )-(4-15 ), for the higher harmonic contributions pro- 
portional to The appropriate set of bulk equations reads, if we use again the 
representation with a scalar potential and a vector potential. 



-df 



afa,(p.(%f-2/.26«4^))] (4-50) 



The first equation determines the pressure contribution Since the pressure 

appears only in the normal stress boundary condition, this is dealt with in the sub- 



sequent section. Next we construct a particular inhomogeneous solution of Eq. (4-52) 
for the vector potential The most general ansatz necessary reads 

^(2,2)inhom ^ ^ 9, (l^ivt?^(z) ^mC,M + •^iv°]j?.'^(^) CiV^M + C.C. 

N,M ij 



(4-53) 

Here, summation over all relevant modes i,j is implied (e.g. G {1,2,3}, 

G and i = j = 1 for hexagons, squares, and rolls, respectively, Figj2]) 



as well as over right and left traveling waves {N,M} G {R,L}, cf. Eq. (2-15). Sub- 



stituting this ansatz into the dynamic equations and matching the coefficients with 



the inhomogeneous contributions of the vorticity equation (4-51) yields the func- 
tions ^NA^iji^) ^^'^ '^NMiji'^)- Since their general form is extremely bulky, we list 
in App|c|only 'l^'^M^iz) and '^'j^Miji^) hexagonal {ij = ji = 13 = 23 = 31) and 
square patterns {ij = ji = 15) as well as for stripe solutions {ij = 11). 

The general solution is the sum of the particular inhomogeneous and a general 
homogeneous solution, •Z'f = !Z^(2'2)-hom ^ ^(2,2)hom_ ^^^.^^^ effective 



tangential boundary conditions (cf. App B.2) 



dy ^ ^{F'NMijiiNijM + F'NMijiiNCjm + CC.) (4-54) 

N,M i,j 

with a suitably abbreviated amplitudes F'j^j^j-. The special form of the right hand 



side is obtained, if in Eq. (B 5) the first order expressions for the variables are 



explicitly put in. Substituting the inhomogeneous solutions t^?^2,2)inhom .^^^ ^^^^^ 
a modified boundary condition for the homogeneous solution results 
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dy ^ "^^{FNMijiiNijM + FNMijiiNCjm + C.C.) , (4-55) 

N,M i,j 

since the inhomogeneous solution does not satisfy the boundary condition. In par- 
ticular, on the right hand side the inhomogeneous part of the boundary conditions 
at z = is modified 



FnM' 



1-3 



with 



Tiinhom 
NMij 



jhl 



f?l I riinhom 

^ NMij "T ^ NMij 



CiN^. 



.jM 



(4-56) 
(4-57) 



Similarly one obtains the y component of the tangential boundary condition starting 
from Eq. (B 6). 



Now the general solution of^P^'"^^^""^ , Eq. (4-55), can be obtained using an ansatz 
similar to that used for the solution of the magnetic potential, Eq. (4-8) 

^{2,2) ^ Yji'fNMij e^^^^^kciiN^jM + 0NMij e^^^'^kcCiN^jM + c.c.) (4-58) 

N,M i,j 



(2,2)hom 



N,M i,j 



(4-59) 



where, again, the first summation is over right and left traveling waves and the 
second one over the fundamental modes involved. The inverse decay length for the 
rotational flow contributions, q, depends on the angle between the i-th and the j-th 
mode 



Qiij 



^lij + 



(0)i2 



fl2 + V2 



(0) 



(4-60) 



and accordingly q2ij by substituting /cgjj for k]^^j in Eq. (4-60), where kuj is defined in 

Eq. (4-5). Here, is an abbreviation for the Fourier transformed time derivative 
and takes the values ioj^'^'^ + a^^\ a^^\ and — iw*^^) + a^^") when applied to RR, RL or 
LR, and LL modes, respectively. The bulk equations and boundary conditions are 
fulfilled for the amplitudes 



hom 



{FRRij - '^iJ'2hijkc^RRij) 



and 



"^RRijiiRijR 



llii + ^ 



lij 



kcDf\iR^jR 



(4-61) 



(4-62) 



Frr 



iiRijR - 2e^'^5,4^) 



llij + k^ 



lij 
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For the last expression we explicitly used the kinematic boundary condition for the 
second perturbative order, Eq. (4-19). In AppjC] these solutions for the flow poten- 
tials are specified for hexagons, Eqs. ( [CfT2l ) and ( [CfTsl ), and squares, Eqs. (|Cll3l) 
and (C 16). The amplitudes with a tilde are obtained from those without one by 



replacing kuj or quj by k2ij or q2ij, respectively. For (fRRij^iRCjR this leads to a 



denominator ~ /c2ij, which vanishes for i = j according to Eq. (4-6). Nevertheless 



all physical quantities derived from that potential, like velocities and strain com- 
ponents, stay finite. The amplitudes in Eqs. (4-61) and (4-62) for the RL and LL 



(instead of RR) components are obtained by choosing the appropriate expressions 
for qiij and d[^\ according to the rules given above. The only remaining condition 
not yet satisfied is the normal stress boundary condition, which we will discuss in 
the next section. 

4.3. The normal stress boundary condition 

To find the solutions to the hydrodynamic bulk equations ( 4-13||4-15 ) , it was not 
necessary to use the normal stress boundary condition. The same situation appears 
in the derivation of the linear eigenvectors. There, substituting the eigenvectors into 
the normals stress boundary condition yields the dispersion relation restricting the 
linear solution to those with a specific uj{k) relation. The second order normal stress 
boundary condition, as will be shown below, leads to the determination of M^^\ 
the first correction to the control parameter entering the final amplitude equation in 
linear order. 



The second order norm al stress boundary condition has been derived in App. B.2 
and is gi ven as Eq. (B 3). It consists of two parts, one is proportional to 



Eq. (4-63 ) and the other to ^^^^ . The latter equation can easily be fulfilled by splitting 
the pressure = ^(2,2)-B _^p{2,2)s ^^^^ part, ^(^'2)^, that is determined by the 
bulk equation Eq. (4-50) and the other, p^^'^)*^, by the ^^'^^-boundary condition. This 
ansatz works, if Ap^"^^ = in the bulk. Indeed, ^(^,2)5 ^ ^.^.^kujz ^ ^,^*^k2ijz 
leads to the required result. This additional pressure contribution is due to the 
inhomogeneities arising in the normal stress boundary condition, in particular the 
one due to surface tension. Since the surface tension always acts normally to the 
surface, this is the only point, where it can enter the nonlinear dynamics. It just 
contributes to the Laplace pressure, which is proportional to the curvature of the 
surface, a quite intuitive result. 

However, this additional pressure contribution is of no importance because of two 
reasons. First, the pressure always enters linearly the hydrodynamic bulk equations 
and therefore it will never give rise to inhomogeneous contributions, which have 
to be accounted for by Fredholm's theorem. Second, the pressure enters only the 
normal stress boundary condition, which is actually the governing equation for the 
appropriate pressure contribution in the next order. In addition, also p^^'^)^ is not 
needed in the following and we refrain from showing it here. 

The situation is different for the first part of the normal boundary condition 
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= J^M(i)M,e(') . (4-63) 
1+^ ^ ^ 

It serves to determine the yet unknown control parameter M^^\ which defines the 
expansion parameter e, on which the amphtude equation concept is based on. In 
contrast to bulk instabilities, where M^^^ follows directly from Fredholm's alterna- 
tive, here we have to employ the normal boundary condition, since the Rosensweig 
instability basically is a surface instability. The same is true for the Marangoni 
instability, where again the driving force of the instability is not contained in the 
bulk equations, but acts purely at the surface. In some previous discussions this 
problem was circumvented by using a scalar product artificially implementing the 
driving force into Fredholm's theorem. This special scalar product made use of the 
fact, that the free boundary was treated as undeformable. In the presence of a 
deformable surface, however, this specific scalar product seems to fail. 



Satisfying the normal stress boundary condition (4-63) provides us with the 
necessary relation between the control parameter M^^-* and the scaled gro wth rate 
a^^\ The actual calculations to get these quantities are displayed in App 
here the results are given 



D.l 



and 



a 



(1) 



i^2(l + /U) 



and Lo 



(1) 



0. 



(4-64) 



The fact that uj^^^ vanishes states, that the instability remains stationary and ex- 
cludes possible soft mode oscillatory branches beyond the linear threshold. For the 
slow growth rate a^^^ we obtain the physical result that the growth is the faster the 
farther one is beyond the linear threshold and it is the slower the more viscous the 
medium under consideration is. We also observe that the elastic contributions in 



Eq. (4-63) cancel upon substituting the solutions of the eigenvectors. This, on the 
one hand, states that Eq. (4-64) applies to ferrofiuids and ferrogels, alike, and on 



the other hand it states that the growth process is solely given by the dissipative 



mechanisms in the system under consideration. Eq (4-64) additionally tells us, that 



the boundary behaves qualitatively di fferen t with respect to the temporal properties 
when compared to the bulk (cf. Eq. (4-28)), since it does not scale with a^^\ This 
qualitative difference is manifest in the kinematic boundary condition which always 
connects the velocity field to the temporal change of the amplitude, as we discussed 
already. It is therefore reasonable to compare the scaled time derivative from the 



bulk with the time derivative in Eq. (4-64) 



In order to com bine t he results from the surface with the solvability condition of 
the bulk equations (4-28), we need to rewrite the growth rate 

in dimensionless 

form. By multiplying Eq. (4-64) with the typical (linear) time scale tq = V2kc{pG + 



A*2^c) ^ and by defining tqit^^) as we obtain 



2{l+ fi){pG + fi2kc 



(4-65) 



and by adding the solvability condition from the bulk equation (4-29) with the one 



from the surface (4-65), we finally end up with a rudimentary form of an amplitude 
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equation for the second order 

In the last step we expHcitly assumed that the dimensionless time derivatives at the 
surface and in the bulk are of the same order. By adding the two subsystems we 
therefore accounted for the singular behavior of the kinematic boundary condition. 
By now we have solved the second order problem completely, with the amplitudes 



of the critical modes satisfying Eq. (4-66). 



§5. Third order 

With the complete solution of the second order problem at hand we can now 
discuss the third order, in order to obtain the desired amplitude equation. As in the 
second order, the solvability condition consists of two parts. One due to Fredholm's 
theorem and one that guarantees the normal stress to be compensated at the bound- 
ary. However, we will have to find solutions of the third order problem proportional 
to the main characteristic modes, only. 

5.1. Magnetic contributions proportional to ^^^^^ 

We again start with the magnetic contributions and restrict our discussion to the 
parts proportional to the main characteristic mode The differential equations 
for the scalar potentials of the distortions to the magnetic fields read 

Z\^(3.1) = and Z\^(3,l)vac ^ Q (5.^) 

with the corresponding boundary conditions (at z = 0) given by 

g^^(3,l)vac _ 9^^(3,1) = M(2)ay^(l) (5-2) 
g^^(3,l)vac _ 5^^(3,1) ^ (5-3) 
^^^(3,l)vac _ ^5^^(3,1) ^ . (5-4) 

The solutions of this set of equations is obtained following the lines of the second 



order calculations, Eqs. (4-9) and (4-10), leading to 



^(3,1) = and ^(3,l)vac ^ ^ (5.5) 

The contributions due to the higher harmonic modes could in principle be cal- 
culated in the same way as in the second order. However, these contributions again 
contribute only to the pressure offset and are therefore of no importance for the 
amplitude equation. 

5.2. Hydrodynamic contributions proportional to ^^^^ 

The complete set of hydrodynamic bulk equations for the hydrodynamic vari- 
ables reads in third perturbative order 

^ + _ 2^25,eg^ - U2 {djdivf + djd.vf) (5-6) 
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a(0) i3) 



r,{2) (1) r,{l) (2) 

-dAp^"^ +P^^ -2,.^^ -2,,^^) 



o(2) (1) ;^(1) (2) (1)^ (2) (2)^ (1) 



(3) 



0. 



(5-7) 
(5-8) 



We restrict our attention now to the contributions proportional to The 
hydrodynamic equations, written in terms of the flow potentials, then reduce to 

A^^''^ = (5-9) 
pZ^a(0)^(3) + ^p(3) ^ _a(i)p/^^(2,i) _ a(2)pz^^(i) (5.10) 

-P2Ad?^^^^)-p[df^M'^^^). (5.11) 

To find the solutions, we follow the same lines as in the previous order. The 
particular inhomogeneous solutions for the vector potential read 



l^(3,l)inhom 



^zba 



P2+P2 
QP2 



d 



V/i2 



ze 



qz 



for {a,b} £ {x,y} while the inhomogeneities in (5-10) are compensated by 

^(3,l)inhom ^ _paW^(2,l) _ pd^^^ ip^ . 

The general homogeneous solutions take the form 

^|3,l)hom _ _5^i^(3,l)g'?z^{l) ^ ^^3,l)hom _ ^^i^(3,l) g^^^Cl) 

and (^(3,1) = ^(3,i)gfc..^(i) ^ 
where the amplitude for the vector potential is given by 



(5-12) 
(5-13) 



The homogeneous solution for the pressure reads 



(3,1) 



P2 + P2 



(2) 



(l)l2^ 



(0) 



P 



(3,l)hom ^ _pqW ^(3,1) 



(5-14) 



(5-15) 



(5-16) 



and upon exploiting the kinematic boundary condition we obtain the amplitude for 
the scalar potential 



(3,1) 



(2) , [dl 



(l)l2^ 



>2(g^ + ^c) 



or + 



(0) 



(5-17) 
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As in the second order, the normal stress boundary condition is not used when 
deriving the solutions. Again, it allows to calculate the linear contributions to the 



amplitude equations. This is done in App. D.2| with the result 

li2kc r5(i)i2.m A;e/x(2M(2)M, + [M(i)]2) 



+ 



pG + ji2K 



2{l+ p){pG + ^i2kc 



(5-18) 



where dj^ = roo"^^-*. Note that the second term is absent in a ferrofluid without 
elasticity. 



§6. Amplitude equation 



We are finally left with satisfying Fredholm's theorem for the third order bulk 



hydrodynamic equations. The general solvability condition for the equations (5-6) 



(5-8) reads 



{v, I p d?vf^) + (e., I r^eljO + {v. I Pdrvr>) + (e., | 



,(2) fl) 



W,,(2,l)\ 



^{1)J2,1)> 



(6-1) 



+2p2{vi\dj{e 



(1) (2,1) , (2,1) 



(1)^,(2,2) ^^(2,2)^,(1) 



(1) (2,2) , ^(2,2)^(1) 
jk ki 



(2,2) 



,(2-2); 



where we already separated the contributions from the second order eigenvectors 
that are proportional to ^(^^ from those proportional to ^^'^K The first two contri- 



butions on the left hand side of Eq. (6-1) can be discussed in the same way a s the 
equivalent terms in the second perturbative order by replacing in Eq. (4-23) a^^^ 
and uj^^^ with (j^^-* and w^^-*, respectively. Thus these contributions yield the scaled 
dimensionless time derivative d!^^ = a^"^^ /a^^^ for the bulk part. The third and the 



fourth contribution on the left hand side of Eq. (6T) can in principle contribute to 



the second time derivative, since the second order eigenvectors w^'^^ and e^j'^^ are 



proportional to d'^'^ (cf. Eqs. (4-43)-(4-49)). Discussing the last contribution first 



we obtain upon exploiting the result of the second order, oj 

N 



i=l 



(1) = 0, 
(0)l5^ 



(6-2) 



This contribution is at least of the order [fjC^^J^ and therefore vanishes in the limit 
of a static instability. Similarly, the contribution due to {vi \ d'^^^^'^'''^ 



is at least of 



the order [cr^'^^]'^ and can also be neglected. Let us now focus on the right hand side 



of Eq. (6-1) and discuss those contributions first that are due to the eigenvectors 

(2,1) 



e ' ' and v. 



(2,1) 



of the second perturbative order, which are proportional to the main 
characteristic modes These contributions involve the combinations of three 
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amplitudes ^^^^ and due to the lateral integ ration they therefore remain finite only 
in the case of hexagons. If we use Eq. (4-28) to substitute e.g. Ci^20'^^^S,3 by —kccr^'^^ 
I ^1 Pi we finally obtain 



(2,1)\ 



64 



(6-3) 



Note, that this term only contributes to the cubic coefficient for the hexagonal pat- 



tern and vanishes for any other pattern. All other contributions in (6-1) involving 



the eigenvectors ' or e^^ ' are at least of the order [cr*-'^^]^ and vanish in the static 
limit. The remaining contributions involve the eigenvectors of the second perturba- 
tive order that are proportional to the higher harmonics Since their analytical 

expressions are bulky, the corresponding contributions to the cubic coefficients have 

(1) (2 2) 

been calculated with Mathematica. For the term {iij \ dl e\,' ) one has to exploit 



Eq. (4-28) in the same manner as done for Eq. (6-3). The final results for the cubic 



coefficients A' and B'{ 
consideration, by 



) are given, for the different regular surface patterns under 



B' 



A' 
27r/3) 
= vr/2) 



184/X2/cj? 

(1256315969/10368 - 69828^3) /i2A:e 
(31831/2 - 11072\/2)^2fcc 



(6-4) 
(6-5) 
(6-6) 



and the solvability condition in the third perturbative order which is due to the bulk 
equations can be written for the hexagonal pattern as 



A' 



16/Li2^c 



111 I'll 



(6-7) 



with all its cyclic permutations 1—7-2 
case of the square pattern 



B'{9,,=2n/3) 2 , |p |2w 
32^2k. ('^'' +'^^' 

-)■ 3 — 7- 1. Correspondingly one finds in the 



A' 



16/U2/Cc 



- 2P B'{e,,=7r/2) - 

?i 4i 57^ — ; ?5 4i • 



(6-^ 



From those equations (6-4)-(6-8) it becomes clear that the dependence of the cu- 



bic coefficients on the material parameters is solely given by the characteristic wave 
vector kc- Thus they are independent of the elastic shear modulus and the mag- 
netic susceptibility. The same is true for the quadratic coefficient as observed in 
Eq. (4-29). This behavior could have been anticipated by inspecting the general 



expressions for Fredholm's theorem (Eqs. (4-20) and (6-1)). The lowest order in the 



expansion with respect to d^^^ is always proportional to the shear modulus fi2 (since 



the adjoint strain field, Eqs. ( 3-8||3-13 ), is proportional to the shear modulus) which 
therefore cancels in Eqs. (4-23) and (4-24). This behavior is due to the assumption 



of linear elasticity. Similarly the assumption of a linearly magnetizable medium and 
neglecting magnetostrictive effects results in cubic coefficients that are independent 
of the magnetic susceptibility. 
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Adding Fredholm's theorem in the third order expansion (6-7) to the correspond- 
ing solvabihty condition from the normal stress at the boundary (5-18), we obtain 
for the hexagonal pattern 



2{pG+H2kc 



_ A:c/x(2M(2)Mc+[MW]2) . A' 



4il+fi)ipG+p2kc) 32p2kc 

B'{eij = 27r/3) ^ ,2 , |P |2^P 



111 I'll 



(6- 



where we assume, as done in the second order, that the scaled time derivatives at the 
surface and in the bulk are the same, because of the kinematic boundary condition. 

Recall now the results for the hexagonal pattern that we obtained from the 
solvability condition in the second order, Eq. (4-66) 



2{1 + n){pG+p2kc 



I^C U A* 



-?2?3 • 



(6-10) 



If we follow the standard methods and multiply the third order equation (6-9) by 
and the second order equation (6-10) by e^, we obtain 



IJ-2kc 



A:c/i(2e2McM(2) + e2[M(i)]2 + 2eM^M^^^ 



A{1 + p){pG + ^l2kc 



<1 



2 P* 



A! 



32fi2k, 



^3, ,,2, i?'(% = 2vr/3^|^.^|,^|^.^|,^^.^ 



■e l?i| C 



64/i2fcc 



(6-11) 



By the series expansion of the magnetization, Eq. (2-12), we can write 

2 



- = [Mc + eM(^) + e^M^^) + . . . 



(6-12) 



and define the control parameter e in the usual way as the relative quadratic deviation 
from the critical value (of the magnetization) 



(6-13) 



Substituting the series expansion of the time derivative in terms of e (cf. Eq. 2-16) 
and using the standard scaling 

ekc^/Aii , (6-16) 



65« + e'dP 



(6-14) 
(6-15) 
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Ah, A 




Fig. 3. Bifurcation diagram: The stationary amplitudes of the hexagonal pattern, Ah (the 
in Eq. (6-17l) and the square pattern, As (the in Eq. \6-19 \) as functions of the control 
parameter; the special values of the latter are explained in the text. The topology of this 
diagram is similar to that obtained by the energy methodEi) 



the amplitude equation can be written as 



1 



9t6 + 7,^X^1 = 



1 



2VA 



4243" 



16 1' 6 



B 



120 



(I6|' + k3|')6, (6-17) 



where we introduce the dimensionless parameter 6 = fi2kc{pG + ^2^0) ^ and where 
the abbreviations A and i?i2o are given by 

A' _ , „ B'( 



A 



5.750 and B 



120 



:27r/3) 



64/i2A;3 



3.544. 



(6-18) 



Starting from Eq. (|6-8|) instead of Eq. (|6-7|) we obtain the corresponding amph- 

B90 



tude equation for the square pattern 



161' 6 



2-.. — - A 
where the cubic coefficient Bqq is analogously given as 

B'(% = 7r/2) 



90 



64/^2^:3 



4.021 . 



(6-19) 



(6-20) 



The fact that the linear contribution on the right hand side of Eqs. (6-17) and (6-19) 



is only proportional to the control parameter e justifies a posteriori our choice of the 
typical time scale tq. 



Let us first consider the static solutions of Eq. (6-17) summarized in Fig. 3 The 



quadratic contribution gives rise to a transcritical bifurcation from the flat surface to 
a hexagonal pattern at the linear threshold!^ A bistable regime exists for negative 
control parameter values e with its lower boundary given by 
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The solution for the hexagonal pattern takes the form 
where the magnitude of the amplitudes reads 



Cile^"^' fori e {1,2,3}, 



1 + yJl + S{A + 2Bi2oYe 



4V^(1 + 2Si2oM) 



(6-22) 



and where the phases have to fulfill the condition = 0. 

Investigating the values of the cubic coefficients we realize, that B120/A < 1 
indicating that the hexagon solution is always stable with respect to stripe solutions 
at the linear threshold. Stripes and squares are mutually exclusive pattern and since 
-B90 + 2i?3o < A + 2i?i2o and Bi^q/A < 1^ the hexagons are losing stability with 
respect to squares at the critical control parameter given by 



90 



2B' 



30 



2{A + 2^120 - ^90 - 2B: 



30) 



(6-23) 



where the cubic coefficient B^q ~ 4.188 describes the nonlinear interaction between 
the hexagonal and the square pattern. 

The square pattern is stable for control parameters larger than 



A + B, 



£5 



90 



2{A + Bgo - B 



120 



(6-24) 



Since £5 < 63, also a bistable regime between the hexagons and squares exists. 

Let us now focus on the dynamical behavior of the patterns beyond the linear 



threshold. We assume that the hexagonal pattern with the amplitude | | , Eq. ( 6-22 ) 



has developed and disturb it homogeneously in space by a small excess amplitude r, 
I I + The linearized amplitude equation (6-17) for the disturbances r then 

reads 



1 



1 



A 



3 1 + 2- 



B 



120 



A 



r . 



(6-25) 



Substituting the solution ( |6-22 ) in the right hand side of Eq. (6-25 ) it can be simplified 
to — (e/2+ /(2\/]4))r, which is always negative above the linear threshold. This 
reflects the fact that the exponential growth of the infinitesimal disturbances of 
the flat surface above the linear threshold gets nonlinearly saturated by the cubic 



coefficients and a stable pattern develops. Eq. (6-25) therefore takes the form of a 
damped harmonic oscillator which can be solved by using the ansatz r =| r | e^"^ 
with the eigenvalues 



^1/2 



1 



A6 



(6-26) 



where the eigenfrequency i? of the oscillator is given by i? = {e\/ A+ \ )/\/A6. 
These last results are still in dimensionless units. If we choose the time scale 



1^2 1 fJ'2 to compare dissipative and oscillatory processes as suggested by Eq. ( D 9 ) 
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< 




2 3 

dimensionless time T 

Fig. 4. Qualitative time dependent behavior (not to scale) of the surface spikes according to 



Eq. (6-25 1. The time T as well as the amplitudes | \ are dimensionless variables. If the control 



parameter e is slightly beyond the critical threshold the plot can be viewed as the qualitative 
dynamics from the flat surface | | = to the spiked surface | |= 1. 



the eigenvalues read 



Al/2 



V ^^1 VAi^I 



This result is intuitive, since the damping rate is inversely proportional to the dissi- 
pative processes, given by 1^2, whereas the eigenfrequency increases with increasing 
shear modulus. We also realize that the relaxation towards the equilibrium pattern 
becomes faster in a stronger gravitational field as well as for larger surface tensions 
and elastic higher shear moduli of the medium. 

The bifurcation from the flat surface towards hexagons is transcritical and there- 
fore involves a non continuous transition. If the control parameter is slightly above 
its critical value, the still flat surface (at T = 0) can be interpreted as a disturbance 



to the stable stationary solution (6-22). The dynamics towards hexagons from the 
flat surface is then described by equation (6-25) giving rise to an overshoot and a 
damped oscillation towards the equilibrium value (cf. Fig. [4]). 

§7. Usual Ferrofiuids 



In the case of ferrofiuids, the dynamic equation for the strain field (2-2 ) is absent 
and we only retain the continuity equation (2-3) and the Navier-Stokes equation (2-1 ). 
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Furthermore, all contributions to the stress tensor (2-6) that are proportional to the 
elastic shear modulus /X2 drop out. The expansion to the nonlinear regime follows the 
same lines as for magnetic gels. The solvability for the bulk hydrodynamic equations 
for ferrofluids then reads 



{v,\d['\pvf^)) = -{v,\d,{pv^\^^)) . 



(7-1) 



Following the same lines as in ^4.2.1 we obtain for the left hand side of Eq. (7-1 ), if 



we retain the lowest order of the expansion in terms of ^^ 



(0) 



{v^\d\'\pv^^)) 



12p 



N 



(0)1 2 



(0)12 



(7-2) 



Similarly the right hand side of equation (7-1) reads 

-24pa(0)([^W]2 



{V:\d,{pvf\f^)) 



(7-3) 

We realize, that the lowest order in the expansion of Eq. (7-1) in terms of d'f^ is 
at least proportional to [9^*^^]^ due to the deformable surface: The velocity in the 
original and the adjoint space have to be proportional to the time derivative of the 
surface deflection and, consequently, the product {vi \ dj{pv^^\j^^)) is proportional 

to This is also the reason why the contributions to the amplitude equation in 

the case of ferrofluids do n ot co ntribute in the case of magnetic gels. The common 
factor - [cr(0)]2) ^j.^^ cancels and, finally, we end up with 



a 



(1) 



and 



|Cl| =161 =1^3 1 



ai^r' 3 

and the corresponding conditions for all cyclic permutations 1 



(7-4) 



1. As for 



magnetic gels, Eq. (7-4) only exists for the hexagonal pattern, whereas for any other 
surface pattern the amplitudes show no nonlinear interaction in the second order. 

The solutions for the second order eigenvectors can be taken from the discussion 
of magnetic gels in §^4.2.2 and 4.2.3 by simply substituting /U2 = and are therefore 
not shown here. With the solutions of the second order, the third order Fredholm's 
theorem can be fulfilled. The latter reads in the case of ferrofluids 



(^i I pdi 



(2)„,(1)\ 



V 



^(1)„,(2.1)\ 



{vi\pdl 'v^ 
-p{vi\dj{vl''v^ 



(1)^(2,2)^ 



I ^ . (1) (2,1) , (2,1) (1) 

p{vi\dj{vl 'vj +vl 'vj 



)) 



(1)„,(2:2) ^^(2,2)^a)^^ ^ 

,(2,2) 



(7-5) 



Since the analytical expressions for the eigenvectors ' are bulky, the explicit 
calculation of the cubic coefficients has been performed with Mathematica and the 
results are shown below . W e should mention, however, that also in the third order the 
right hand side of Eq. (7-5) is proportional to [Sj*^^]^ and the global factor ([o;^'^^]^ — 
[cj*^*^)]^) can be canceled. 
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The discussion of the normal stress boundary condition in the case of ferrofluids 



can be taken from §^ 4.3 and 5.2 In the second order, the additional condition to the 



amplitudes (4-65) is valid for ferrogels and ferrofluids, alike, and in the corresponding 
third order condition (5-18) we have to substitute /i2 — with the consequence that 
there is no second order time derivative. The typical time scale in the case of 
ferrofluids is then given by tq = f2^c/(pG') which is in accordance with previous 
theoretical discussions EH' The final amplitude equation is derived in the same way 
as in ^for magnetic gels and finally results for the hexagonal pattern in 



5r6 



1 



r?2?3" 



2 3VA^ 
with the control parameter e-^' defined by 



aP6 



B 



120 



fl 



where 



G is the linear threshold for ferrofluids 



(7-6) 



(7-7) 



For the square pattern the quadratic coefficient is absent and we obtain 



5t6 = e"6- 161 
where the cubic coefficients are given by 



Rfl 

-^90 

A^ 



B 



120 



B, 



90 



8.625 
3.150 
4.266. 



(7-8) 



(7-9) 
(7-10) 
(7-11) 



The discussion for the different stable patterns follows the same lines as in ^ 
At the linear onset we find hexagons to be the preferred pattern, which remains 
subcritically stable for control parameters larger than 



9(Afl + 25^20) 



(7-12) 



Since Sgg 



+ 2B^ 
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2Bf2Q and B^q/A^ < 1, where the cubic coefficient ac- 



counting for the nonlinear interaction between hexagons and squares is given by 
B^Q ~ 4.545, the hexagon pattern transforms into a square pattern for control pa- 
rameters larger than 



2(5^0 + 2i?fo) 



(7-13) 



9(^fl + 2Bf,, - Bl, - 2B%Y ■ 
The square pattern in turn becomes unstable again for control parameters lower than 

es = 'i^' + O (7.14) 

' HA^ + Bl,-BU-Bl,y ^ ^ 

In contrast to the nonlinear discussions based on the energy methodll''ll2l' we find 
stable solutions beyond the linear threshold. In particular, our stability boundaries. 
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except the linear stability, do not depend on the magnetic susceptibility. The reason 
for this is twofold. First we neglected magnetostriction resulting in a decoupling 
of the magnetic and hydrodynamic degrees of freedom in the bulk. And second, 
we assumed a linear magnetization law for the superparamagnetic medium. The 
bifurcation scenario in realistic ferrofluids, however, depends on the magnetic sus- 
ceptibility which we trace back to the nonlinear magnetizability of ferrofluids and 
to the fact that the macroscopic parameters, as the viscosity, depend on the mag- 
netization. To capture these effects, one therefore has to start from a macroscopic 
model for magnetic fluids with nonlinear material properties. 

§8. Discussion 

In this paper we have derived the amplitude equation for the Rosensweig insta- 
bility in isotropic magnetic gels based on the fundamental hydrodynamic equations. 
An important step was to find the adjoint linear system of equations together with 
its corresponding boundary conditions in the presence of a deformable surface. Two 
assumptions turned out to be crucial in order to find the adjoint system. Besides 
the dynamic treatment of the Rosensweig instability, the medium has to be consid- 
ered compressible for the adjoining process. The reason for the latter assumption is 
to maintain the symmetry of the stress tensor during the adjoining process. While 
we can assume an incompressible medium after the adjoining process, the dynamic 
treatment of the system of equations turns out to be also important in the discussion 
of the higher pcrturbativc orders. 

With the help of the adjoint system we were able to satisfy Fredholm's theorem 
and to perform a weakly nonlinear analysis. However, due to the decoupling of the 
magnetic bulk equations from the hydrodynamic ones, Fredholm's theorem does not 
contain the control parameter, which enters the boundary conditions, only. 

We solved this problem by observing that the normal stress boundary condition 
consists of two parts. One is proportional to the higher harmonics of the character- 
istic wavelength and merely increases the hydrostatic pressure in the medium. The 
other one is proportional to the main characteristic wave vector and serves as an 
additional solvability condition providing the dependence between the scaled growth 
rate and the control parameter. Both solvability conditions show qualitatively dif- 
ferent behavior in the static limit. While the solvability condition obtained from the 
normal stress boundary remains finite, the bulk contributions scale with the linear 
growth rate. The latter behavior is mediated by the kinematic boundary condition 
and has been taken into account while combining both solvability conditions into 
one. Furthermore it reveals the fact that both states, the initial flat surface and the 
final spiked one, are motionless states where the velocity field vanishes identically. 
This does not mean, however, that the Rosensweig instability can be simply treated 
as a static, energetic one, and the present work has shown why. 

While combining the bulk solvability condition with the normal stress boundary 
one has some freedom to choose the relative weight of the boundary with respect 
to the bulk via the two differently scaled time derivatives. It seems most natural to 
weigh these single contributions equally with respect to each other. We note that 



Amplitude Equation for the Rosensweig Instability 



31 



this has also been done imphcitly, for example, in the nonlinear discussions using an 
extended scalar product!^ 

Upon combining both parts of the second and the third order solvability condi- 
tion, then following the standard procedure, we obtained a set of amplitude equations 
for the special cases of stripes, squares and hexagons. The latter contains a quadratic 
coefficient that renders the bifurcation from the flat surface to the hexagonal pattern 
transcritical. The calculated cubic coefficients additionally reveal that at the linear 
onset hexagons are the stable surface pattern. For high magnetic field strengths 
instead, hexagons become unstable and a square pattern develops. Both transitions, 
from the fiat surface to hexagons and from hexagons to squares, involve bistable 
regions. We obtain qualitatively the same results in the case of ferrofiuids, where 
the derivation of the corresponding amplitude equation and the determination of the 
nonlinear coefficients has been discussed in ^ We note, however, a qualitative differ- 
ence regarding the temporal derivatives. For magnetic fiuids the amplitude equation 
is first order in time, while for magnetic gels it also contains a second time derivative 
reflecting the elasticity of the gel. Amplitude equations that are second order in 
time are well established in problems related to the buckling of plates and shells. 
For example, Lange and Newell^*^' have analyzed in detail the post-buckling problem 
for thin elastic shells. Along the same lines one has already derived phase equa- 
tions, the analog of hydrodynamic equations for large aspect-ratio pattern-forming 
systems, containing first and second order time derivatives for elastic systems under 
the influence of an external 

The results for the static patterns in this article are in qualitative agreement 
with the bifurcation scenario obtained with the energy methodIl3 The cubic coef- 
flcients in the present paper, however, are independent of the elastic shear modulus 
and the magnetic susceptibility. This is due to the assumptions of ^ where we 
modeled the magnetic gel as a linear elastic and a linearly magnetizable medium 
and where we neglected magnetostrictive effects. Within the energy method, where 
the same approximations have been used, the fourth order coefficients (these coef- 
ficients qualitatively correspond to the cubic coefficients in the present treatment) 
showed an inverse proportionality on the control parameter e. As we analyzed in 
Ref. |11( ), this is due to independently minimizing the energy density with respect 
to the higher harmonics and the main characteristic modes in that method. This 
dependence has been omitted in the subsequent discussions of the energy method 
for simplicity rendering this approach valid in the asymptotic limit of a vanishing 
magnetic susceptibility only. In retrospect this minimization procedure and the sim- 
plification afterwards appears to be unsystematic. The results of the present article, 
however, are valid for a finite magnetic susceptibility and for finite shear moduli. 

In our discussion of the nonlinear properties of the Rosensweig instability we 
assumed spatially homogeneous patterns with no long wavelength variations. By 
additionally rescaling the spatial coordinates in the same way as the time coordinate 
( 2-16 ), one could also implement these possible variations in space. As a consequence, 
the amplitude equation additionally contains derivatives of the amplitudes with re- 
spect to the scaled spatial coordinates. For typical nonlinear differential equations 
these linear contributions to the amplitude equation can be obtained systematically 
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by a standard method exploiting the hnear properties of the systemllSf'I^ll' 

In the case of the Rosensweig instabihty, however, we additionaUy have to take 
into account the deformability of the surface and along with it the kinematic bound- 
ary condition. If the surface deforms, also its normal vector n changes in the course 
of time, which we took into account in our previous discussions by explicitly ex- 
panding the latter in terms of the surface deflection ^. All the different orders of 
n involve gradients of the surface deflection ^ as can be seen in Eqs. ( |A[2||A[4"| ). 
Upon rescaling the spatial coordinates we also must expand the gradients appear- 
ing in n in terms of e, which leads to additional contributions to the higher order 
boundary conditions solely due to the large scale spatial variations of the normal 
vector. These contributions are not contained in the linear dispersion relation and, 
of course, cannot be implemented into it by any means, since the dispersion relation 
only considers the linear properties of the system of equations and therefore assumes 
a still flat surface. One rather has to expand the set of boundary conditions with 
the scaled spatial coordinates from the beginning. The contributions to the second 
spatial derivative in the amplitude equation may then be separated into those due 
to gradients in the stress tensor (for example djVi), which are the ones that follow 
directly from the dispersion relation, and those solely due to the deformability of the 
surface. Furthermore we have to evaluate the boundary conditions at the physical 
boundary, z = ^. In our calculations we accounted for this fact by expanding the 
eigenvectors in term of ^ around z = 0. This again involves gradients with respect 
to z, which have to be rescaled as well and which are not contained in the dispersion 
relation. Additionally, one has to expect contributions to the second order spatial 
derivatives in the amplitude equation that are due to the bulk equations. In the case 
of the scaled time derivative we showed that possible contributions due to the bulk 
scale out in the static limit, but it is not expected that this is also the case for the 
spatial derivatives. In conclusion, it is rather obvious that spatial gradients enter the 
amplitude equation in the form of diffusion terms, although in the present case it is 
very cumbersome to derive the appropriate coefficient and requires a rather lengthy 
new calculation. 

The amplitude equations have been derived using the critical value kc for the 
transverse wavevector. Above threshold, however, a whole range of k values around 
kc are allowed for the patterns (Busse balloon) and the most unstable mode may have 
a transverse wavevector somewhat different from kc- The explicit determination of 
those properties is again rather unwieldy. There is also the possibility of secondary 
instabilities, which however we have not looked at. 

In addition to the static properties of the surface patterns, the analysis in this 
article provides us with nonlinear dynamical processes. We obtain the typical first 
order time derivative that describes the growth of the surface spikes beyond the linear 
threshold but that also accounts for the dissipative processes in the medium. The 
typical time scale of the growth (or relaxation) processes increases for increasing 
viscosities and becomes smaller for increasing surface tension and shear moduli. 
Additionally, however, we find a second order time derivative in the case of magnetic 
gels. 

Throughout our analysis we treated the system dynamically and allowed for the 
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static limit in the very end. The way we discuss the Rosensweig instabihty there- 
fore hardly differs from discussions of oscillatory instabilities. When considering 
amplitude equations for oscillatory instabilities, Coullet et alHH demonstrated that 
depending on the relation between the cubic coefficients either propagating or stand- 
ing patterns are obtained. In our analysis we realized that the imaginary parts of 
the scaled time derivatives, w^^^ and oj^'^\ vanish and we can therefore exclude such 
oscillatory soft mode instabilities above the linear threshold. 

The analysis in this article elucidated the main aspects of the underlying mech- 
anisms that lead to the Rosensweig instability. But it also unraveled that for a 
better quantitative understanding additional phenomena have to be taken into ac- 
count. Two nonlinear properties have been neglected. The nonlinear magnetization 
behavior, that already affects the linear threshold, and nonlinear elastic properties. 
Additionally, the magnetostrictive effect might influence the bifurcation behavior in 
magnetic gels. 
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Appendix A 

Magnetic boundary conditions 



In this appendix we derive the magnetic boundary conditions for the second and 
third order. Since the surface normal n is not constant but depends on the surface 
deflection (as do the distorted field contributions), a higher harmonic coupling to 
previous orders is possible (in contrast to the system of bulk equations). For the 
upcoming calculation it is useful to determine first the fields at the boundary z = ^ 



!z=0 



(2)' 



z=0 



l[i^W^d! + 2^(^)d.)d'p(%=o) (01) 
and accordingly for the magnetic field H™'^ and the magnetic flux densities B and 



-gvac_ fjij^g contributions in (A 1) that are explicitly proportional to ^^^^ or ^^^^ are 
due to the deformable surface. 

As mentioned, the surface normal n, initially directed parallel to the z— axis, 
changes its orientation in course of time as the surface perturbation grows (cf. Fig.[T]). 
To give a proper expansion of the boundary conditions, we additionally have to 
expand the surface normal as a function of the surface deflection ^(x,y,t) 

n = no + e n^^) + e^n^^) + e^n'^^) ([A]2) 
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with the different perturbative contributions given by 

„(i) = -dy^m , n(2) = -dy^i^) (|A}3) 

and n(3)= -dy^i^) - \{d.,^'~')f{dyi'~^)) - \{dyi^^)f (04) 

With the previous considerations on hand, we are able to expand the boundary 
conditions in terms of e. The fact that the normal component of the magnetic flux 
density is continuous at the boundary gives the following condition 

n • (H™" - H) = n • (B^^" - B + M) = n • M ([A|5) 

Consider the linear perturbative order of the last equation 

n(i) . (H™'^ - He) + n(0) • (h(i)™'= - H^^)) = n^^) • M, + n^o) • M^^) ([A}6) 
For the constant contributions (constant with respect to x and y), we find 

while the contributions proportional to n^^^ cancel identically. The corresponding 
expression for the second order contribution to the applied field, j^^^^'^ _ i^'^-* = 

(2) 

Mz , can be obtained straightforwardly. 

The boundary condition for the tangential components of the magnetic field 



(2-9) is given in linear order 

^(1) ^ ^jjvac _ jj^) ^ ^(0) ^ - h(i)) = ([A}8) 



which can be simplified substituting Eq. (A 7) to (with a £ {x,y}) 

/l(l)--/j(l) = -(a,e«)Me (09) 
In the second perturbative order we find 

= n(2) X (H™^ - He) + n(i) x (H^^)™'^ - H^^) - a<p(^)™" + d^'-^^) 

+ n(°) X (H(^)™'^ - 9^(2)vac ^ ^^(l)Q^(l)vac 

-H(2) + 9<?(2) + A:^(i)acp(i)) (|A}10) 

which is simplified in the same manner (by exploiting the results of the previous 
order) to 

+{da^^^^){dz'P^^^"^^ - dz^^^^) - A;e(^)(aa^(i)™= + da'P^^^) = (|A}11) 
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with a E {x,y}. This immediately leads to expression (4-3) used in the main text. 
Finally we deduce for the tangential boundary condition in the third perturbative 
order 



= n(3) X (H™'^ - He) + n(2) x (u^^^'"' - H^^) - 9^(1)™'= + d<pW 

+ n(^) X ^H^^)™'^ - H^^) - 9c^(2)vac _^ g^{2) _^ J^^W Q ^c^(l)vac _^ ^^(1) ' 
n(°) X fn^^)™'^ - H^^) - 9c^(3)vac _^ g^{3) _ ^WQ__g ^^(2)vac _ ^(2)' 

lyfc2^(l)2 _ ^^{2))0^(l)vac ^ (1^2^(1)2 ^ A^^^'V^^'^ ©l^) 



^2 ' " ^ ^2 

where it will be sufficient for our discussion to consider only th e cont ributions pro- 



portional to the main characteristic modes as discussed in {A.2.2 



Along the same lines the boundary condition that guarantees the continuity of 
the normal component of the magnetic flux density is derived. In first perturbative 
order we get 

^(1) . (gvac _ ^ ^(0) . J^B{l)vac _ ^ q (j^^g) 

which is straightforwardly simplified to 

^(l)vac _ ^(1) ^ Q 

For the corresponding condition in the second perturbative order we obtain 

+ n(°) • (b(2)™^ - b(2) - a^(2)vac ^ i^Q^(2) 

+A;e(i)9^(^)™' + ^ik^^^^ d^^^^) (|A|15) 
which is simplified by exploiting the previous order to 

-{^yi^^^){^i^y^^^^ - + + 5^cp(i)) = (|A}16) 

Finally, the third order boundary conditions takes the form 



= n(3) . 


(Br 


+ n(i) • 


(b(2) 


+ n(0) • 


fB(3) 



V2 - - / V2 

where again it will be sufficient for our discussion to focus on the contributions 
proportional to the main characteristic modes 
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Appendix B 

Hydrodynamic boundary conditions 



B.l. Expansion of the boundary conditions 

In this section we discuss tlie expansion of the hydrodynamic boundary condi- 
tions to the second and third order in terms of e. Recall first, that we require the 
tangential stress at the free surface to vanish whereas the normal stress is balanced 
by surface tension (2-7, 2-8). The contributions of the stress tensor to the differ- 
ent perturbative orders are defined by the expansions of the macroscopic variables, 
Eqs. (2-12, 2-13), and by the expansion of the surface normal n, Eq. (A 2). The 
linear eigenvectors of the hydrodynamic set of equations are either proportional to 
gfcz gg^Jii} ^]^g boundary conditions one has to evaluate them at z = ^ and 
therefore an expansion similar to Eq. ( A 1 ) is needed that explicitly accounts for the 
deformability of the surface. The linear order of the boundary conditions is discussed 
extensively in Refs. [4]),20) and is therefore skipped here. 

B.2. The second perturbative order 

In the second order we find as the tangential boundary conditions involving the 
hydrodynamic fields 

= -e(^)a42^2eiy + Md.vP+d,v^^^)] 
+{dy^W)[2^,ei'J + U2{dyvi'^+dJ^'^)] 

-2(9.^(1)) [/.2(eW + 4)) + z.2(5,4') + + pvi%i'^ (§1) 



2/.2e(2) 



yz 



In Eqs. (B 2) and (B 1) the inhomogeneities on the right hand side have been ab- 

(2) — — (2) 

breviated by f^xz and Oyz , respectively. In particular these inhomogeneities are 
proportional to 



The normal stress boundary condition (2-8) in the second order reads 
= -2p2\^f + p[v^^)f - M^B^[{dye)f + {d^e^f] 

-2pH,{dyC^'^){dyCP^^^) + 2/7^ (5,^(1) )(a,c|>(i)™-) 

-2pH,{dxC^'^){dx^^'^) + 2Fr 
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l + /i 



(1) 



1 + /i 



(IB}3) 



Furthermore we obtain for the kinematic boundary condition in second order 

9(0)^(2) + qW^H) + (^(1) . 9)^(1) = „(2) + ^(1)5^41) ([B}4) 

The physical boundary is at z = ^, giving rise to an additional dependence on ^. 



In Eqs. (B 1)-(B 4) such terms have been made explicit (e.g. the last one of (B 4)). 



Thus these boundary conditions are effective ones that have to be taken at 2; = 0. 



What can be realized immediately in the expressions (B 3) and (B 4) is, that 



two qualitatively different contributions arise. On the one hand we obtain the ex- 
pected contributions proportional to the higher harmonic coupling ^(^^-^ of the main 
characteristic mode. On the other hand, there are still contributions proportional to 
the main characteristic mode ^^^^ itself. The latter will allow us to find the linear 
contributions in an amplitude equation even though no explicit control parameter is 
present in the bulk equations. 

To solve the corresponding hydro dyna mic bulk equations, we introduced a scalar 
'P^'^^ and a vector potential ^'^^^ in ^4.2, to discuss potential and rotational flow 
contributions separately. Following the same lines as done in the linear orderjlJ'llSl 
we can translate the boundary conditions into a corresponding set of equations for 
the amplitudes of the second order potentials ^^^^ and We obtain for the 



tangential contributions (B 2) and (B 1) 



usmg = ^yz 
translates into 



at the boundary. The normal stress boundary condition (B 3) 



\ -\- LI 



(2) 



Eqs. ( |B[5HB[7] ) follow from ( |B[2|]B[3] ) by taking the time derivative with respect to 
the fast time scale t^^^ without loss of generality. This is why in Eq. (B 7) only the 
contribution d^^^p^"^^ and no contribution d^^^p^^^ arises, while dj-^^e^^'' gives rise to 

,(2) 



contributions and ~(?) 



(cf. Eq. (4-14)) 



B.3. The third perturbative order 

We take over the procedure of the previous section to the third order. If we use 



the solutions ( 4-36|4-40|4-42 ) of the hydrodynamic bulk equations in second order, 
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the kinematic boundary condition reads 

9(0)^(3) + Qi^)^m ^ ^(2)^(1) ^ (^(1) . ^^^(2) ^ (^(2) . ^)^(l) 

(11^2 

H^'^dA'^ + l^^'^'dU'^ (08) 

The tangential boundary conditions are of the usual structure and can be written as 

2/.2eg) + U2 (9,43) + 0^43)) ^ ^(3) ([B]9) 

2^2eg) + 1^2 {d4^'> + dyvP) = ([B}io) 
where the contributions to the inhomogeneities that are at least proportional to the 

(3) 

higher harmonic couplings are collected in the abbreviation i^lj , similarly as done 
in second order. Taking the time derivative of Eqs. (B 9 B 10) together with (B 8) 
we find 

= 9f )f2g) + 2^2(a«eg) + df^e^) + 2^.2{v^^^^,e^^ + ^eg) fjll) 

- fi2d.dyWr(^^ - fi2{dl - dl)^^^^ + 2/i2a.a,(/?(3) 

= dl'^Qf) + 2/.2(5«4) + d?\^^}) + 2Mvl'^d,e^i] + 9,41^) (012) 



For the normal stress boundary condition we obtain from Eq. (2-8) 



)vac 

~)vac 



= _ ^(2)vac^^^(l)vac ^ ^^(1)5^^(2) _ ^(l)vac^^^(2)^ 

+4^^ + cJTa • n(3) ([Bjl3) 
which by a similar procedure can be written as 

- (2/i2 dy + pGdy)^P + {2fi2d,d, + pGd,)^^^^ + (2/i252 + pGa,)^^^^) 

1 H~ 

^e^i) + + 42) ^^^(1) _ ^(1)5^^(2) _ ^(2)9^41) 

+4') afc6«) + araf^a • n(3) + V^'^ + ^^^g) (014) 

Appendix C 

Eigenvectors in second order 



In this appendix we give the contributions to the eigenvectors in the second 
perturbative order that are proportional to the higher harmonic couplings ^^'^^ . Due 
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to the fact that we have to treat the system dynamicahy throughout ah orders, the 
expressions become tedious and have therefore been calculated with Mathematica. 
In the following the solutions for the hydrodynamic potentials are represented for the 
patterns under consideration, hexagons {9ij = 27r/3), squares {9ij = tt/2) and stripes 
{Oij = 0) as well as for the interaction between hexagons and squares {9ij = vr/G). 

The inhomogeneous contributions to the vector potential, cf. Eq. (4-53), separate 
into a contribution ~ 6^*^'="'"'^)^ and ~e^'^^. For the hexagonal case {ij = ji = 12 = 
23 = 31) we obtain 

^inhom. ) ^ (fcg + g^)(2/i2g^ - 5n2qkc + p[Dl'W^^'^'Di°^ 

^ 3fcgg(/X2fcg + 2^,,q' - p[Dl'^]')e'^wf^ 
{kt-hklq^+Aq''){p2kl-Ap2q' + p[Df^Y) ^ ' 

with {iV,M} G {R,L}. The abbreviation of^ stands for 2ia;(°) + 2cr(°), 2cr(°), and 
-2ia;(°) +2cj(°) for N = M = R, N M , .md N = M = L, respectively. The second 
coefficient 5^]^^™ ^^^^^ 



dQll 



+ 



klq{klp2-2q^P2 + p[D?^?)e''^^D?^ 
(3fc4 - 7^2(^2 + 4g4)(3fe2/i2 - 4g2/i2 + 2p[Df'\^ 

(A:g+g^)(6fc,V2-10fcgg/X2+MV2+29V2+gp[A^°¥)e^^^"^'^'A^°^ 
2(2fc4 - 2klq - 2,klq^ + 2kcq'^ + q^){2klp2 - 2kcqp2 - q^p2+p[Df'^?) 
For the square pattern we get {ij = ji = 15) 

^inhom ^ (fcg + g^)(4^2fcg-10/i2gfc,2 + 2//2g^ + (fcc + g)ppi°¥)e(''+"-)-Ij;°^ 

2(A:4-2g/c3_2fc2g2+2A:cg3+g4)(^2fc2-2/22(?fcc-/i2g2+p[Z)f)]2) 

^ gfc2(2^^g2_^[^(0)]2),2,.^(0) 

4(A;4-3(72A:2+2g4)(2/i2fc2-4/22g2+p[Dj°)]2) ^ 
= l^iv^Jj^^l^) ©4) 
For the stripe geometry, i = j, we obtain 

(A;2 +g2) (6A;2/X2 -4fcegM2 - 2gV2 - pii^f ]2) e('==+9)-i?f 



2(fce-g)(3fc2+4fcc(?+g2)(3A:2/i2-2A;c9/i2-g2/i2 + 2p[Z){°)]2) 



©5) 



and 



inhom 



NMiji^) 
+ Zi 



Aktp2-^\2q^p2-p[Df^?)-kl{2Qq^P2 + MDf^?) j,2^2,.j^m 
4(A;2-g2)2(4^3^2 + ^7/>[A^°¥) ' * 

4(A:,-g)2(A:,+g)3(A:2;i2+2MA2+g'/i2-2p[A^°^]2)j"'e('==+'')^Df) (06) 
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with the numerator Zi being given by 

Zi = {kl + q^) (2Aklqii2-^q^lJi2-2q^p[Df^f + 2Qklq^^i2 + ^klp[DfY 

+Skcq^^i2-^kcqp[Df^^ (07) 

In addition, to describe the interaction between the stripe and the hexagonal 
pattern we need to consider also the case 9ij = 7r/6, {ij = ji = 14 = 36 = 25) 

^mTj = {(A;2+g2)(2g3^2-10A;2g^2+A;,V2(l+3V3)+gp[Z)f ]2 ((C]8) 

+V3MV2 + fec(l-V3)p[Df)]2)}e(^=+'')^Z){°) 

fcgg(fc,V2(2^/3 - 3) - 2(2 - V3){q^p2 - p[Df^Y))e^'^^Df^ 
([2 + Vi\kl - Aq^){kl - q^){[2 + ^^^^2 - ^^P2 + 2p[Df^f) 
where the denominator A'^i is given by 

Ni = 2{kl-q^)[2{2 + ^)ktp2-^{l + V?>)klqfi2 + q^P2-q^p{Df^f 



+ kl[{l + ^)p[Df^f + 2{l-^)q^ft2]+Akcq'^P2-2kcqp[Df^f] (|C}9) 



and 



.Tvinhom 
^ NMij 



with 



^ ((3 + 2V3)fc,V2 + (2 + V3){2q^p2 - pjDiy)) qk^"^ D?^ 
{kl - g2) [( V3 - 2)kl + 4g2] [( V3 - 2)klJi2 + Aq^h - plD^y] 

-N^^ {(A;c+9^) [{3V3-l)klfi2 + Wklqfi2-2q^f^2 

-2qp[Dl'^f + k,V3q'p2-k,[l + V3]p[Diy] } e('==+'')^Df ^W) 

N2 = 2{kl-q^)\2{2-V?,)kt~^2+^{V?,-l)klqp2 + q^P2-q^p[Df^? ([C}!!) 
+2kl [1 + y/i\q^p2 + kl [1 - Vi\p[Df'^f +Ak,q^j:i2 - qkcP[Df^?^ 



For the scalar potential we obtain from Eq. ( |4-62 ) in the geometry of hexagons 
{ij = ji = 12 = 23 = 31) 



SpD)^'q{q+k,){kt-Aq' 

[2klq{l2p2 - 14^2 -3i/2 V 12<7V[A^°¥ 
+2kt (4^2 - 16^29^ + i^2q^Df^ + 2p[Df^f) 
+kl{l^q^p[Df^f + ^q\Vp2+Ap2 + i^2]Df^) 
-Akc {9q'p[Diy- Aq' [3p2 + ^^2!?^^]) 

+fc3(l7g/)[L>{°^]2-96/i29^+52//2g^ + 20g3zy2Z)J°^)} ^12) 
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and for the case of squares {ij = ji = 15) 



fNMij 



AV2pDf^ {kc + q) {kl - 2q^){kl - 2k,q-q^ 



(0)i2 



{0)n 



-Aklq{lQjl2-UiJL2-^y2D 

+2/cc(5gV[A^°¥-V[3/U2 + i^2A^°^]) 
+kl{MDf\^ - 4g2[22A2-14/i2-3i/2A^°^]) 
-kl (35gVpr¥ - V[20/i2 + 15/X2 +3^.21^1°^]) 
-fc,4(l3gp[A(°)]2-4g3[l8A2-22M2-9z.2A^°^])} 

For stripes {i = j) we obtain 



(013) 



{kc-q){kc + q){3kc + q)pD 



(0) 



(12114) 



X {2A:c(2A:c-g)p[I)J°^]^ -z^2(A;c-g)(3fcc+'z)(3A:2-2A:,g+(72)Z)^ 
+ {kc-q){3kc + q) [kl{6fl2-7^l2)-3q^|^2 + 2k^q{fl2 + 2^l2)]} 



(0) 



The homogeneous contributions to the vector potential follow from Eq. (4-61). 
The solutions for hexagons {ij = ji = 12 = 23 = 31) read 



horn 
NMij 



(0)l2^ 



[8kl{k, + q)pDf'\p[Dl^^f-q{2k, + q)){f,2kl-4f,2q'+p[Dl 

\^4:kcqi2kc + q)pl [kt (6/22 - 7/i2 ) + 2klq{fi2 - 4/X2 ) - Sk^q^ {p2 + 2p2 ) 

+A:2g2(i3^2-24/i2)] 
-V2fi2kc{kc-2q)q{2kc + q){kc + 2q){2kl-kcq+2q^)D_ 



(0) 



+/l2p[I?rf [4g^/i2 + 3A:3g2 (^2 - 2/i2) + A^e'g^ (79/i2 + 8/X2)] 

+2kcp[D^^^fp2V2{kl-2kcq-bq^){3kl-k,q+2q^) 
V[A^°¥[5/i2g^ + A:e^(17/i2-8M2) + M'(^2-4/i2) + A:2g(17/i2+4^2) 
+2fcep2^2[I?r¥(3A:,'-M+2g')-p=^[A^°^]^(A:c-g)} 

and for squares {ij = = 15) we obtain 



(015) 



horn 
NMij 



8pDfhl{kc + q) ( (A:,^-2M-g^)/i2-p[I?n 



2\ ~ 



l(0)l2 



xf4/i292-p[Z)f)]2-2^~2fc.^ 
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X { 8A;c/i2 {kl-2k,q-q^) [2k,fL2 {'ik, + q){kl-2q^) 

-fi2{7kt-7ky + 8k,q^-6q^)] 
-%V2hkc{kl-2q^){kl-2kcq-q^){2kl-kaq+q^)Df^ 
+2'^2p[Df^f[kl{'i?,Ji2-2Q^i2)-2q^fi2-2kcq^{bfi2 - 7^2) 

-klq^{2<dil2 + l2^i2)+klq^{2{)pL2-Slfi2) + /c^g(36;U2-15/i2)] 
+Afi2kcPV2{'ikc-bq){k^ + q){2kl-kcq+q^)[Df^f 

V[A^°¥[59^/i2 + A:,^(23/i2-12/U2) + M'(/i2-4/X2) + fccg(ll/i2 + 8/i2)] 
-Ak,p'v2[Df^f{2kl-k,q+q^) + (k,-q)p\Df^f} ([Cjie) 

And for stripes 

^Nm, = - [8klpDf\k, + q) [{k,-q){'ik, + q)J,2+p[Df\^)]~^ 

x\A~^2kc{kc-q){?>kc + q) [fe^(6/i2 - 7 P2)-q^ P2 + 2kcq{ll2 + 2p2)] 
-Afl2kci^2D?\kc-q){3kc + q){3kl-2kcq+q^) 
+p[Dl^^f[fi2i31kl + 3klq + 5k,q^ + q^)-4k,i5kl-4.kcq+q^)] 
-Ak,pu2[D'l''>]H3k^c-'2hqW)+p'[Df^]\k,-q)] ^17) 

With these potentials the components of the vefocity field can be calculated 



straightforwardly. To determine the components of the strain field via Eq. (4-14) 



however, the inhomogeneous contributions —vi^^dke[^^ have to be calculated addi- 



tionally 

vi'^d,e(^^ = (^2^{(1 - + (^c + q')ky'^'' 

-2kl{kl + q^){kl - 2qkcCOs9,, + q'')e^''^^i>]^iNCjM ([CjlS) 



_^(0) a 

^^5,6^ = - cos%) [Aktq'e^'^^ + kl{kl + q')e'''^^] 

-2kl{kl + q^){kl - 2qkccos0^j + q^)e^''^+'>^'}^iN^jM ^W) 
_n(0)|'i,2 ,1.2 \ 

- * f^, _r,^^n.MAl,2^2 2qz . (i.2.2^2 2k,zl 



-2(A:2 + g2)(^2 _ 2gfc,cos% + (z2)e(^+^'=)^}ei7ve,M (020) 

_r)(0)^7.2 I 1,2 \ 

-2{kl + - 2#,cos% + g2)e('?+'==)^}eiivejM (021) 
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1 — cos 9ij ) 



,,(1);^ ,(1) _ iDf'\ki,, + k,,^){e,+q^ 



[2klqe"^' + h{kl + q^)e^''''] 

{kc + q)[2k^,-qkc{l + cos9^j)+q^- kl cos %] e^^^^''^] ^n^jm 
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Q22) 



,(1) 



iDf\h,y + k,,y){klW) 



yz 



1 — cos / 



(lCj23) 



x[2klqe^'i' + kc{kl + q^)e^^^'] 
-(kc + q) [2kl-qkc{l + coseij)+q^-kl cos e^^^+'^^^jeiTvC^M 

where we have displayed only the iiN^jM contributions to —v^k"* dke^^j^ ■ The contri 



butions ~ ^iN^jM "^^^ be derived from Eqs. (C 24) - (C 18) by the replacements 
iiN^jM — > iiN^jM and kj — > -k^, cos% — 



cos Oij 



(IC}24) 



Appendix D 

Satisfying the normal stress boundary condition 



D.l. The contributions in the second order 

In this section we present the explicit derivation of the surface contributions to 
the amplitude equation that are due to the second order surface boundary condi- 
tions, in particular in the case of the Rosensweig instability due to the normal stress 
boundary condition. We restrict ourselves to the contributions in the normal stress 
boundary condition that are proportional to the main modes This is exactly 

the part the gives the relevant condition, while the contributions proportional to .^^^^ 
are compensated by a pressure offset and do not give rise to additional restrictions 

(cf. Q. 

Start ing from the general second ord er no rmal stress boundary condition 
Eq. (B 7), its .^'•^^ part has been given in Eq. (4-63). Using the expression for the lin- 



ear eigenvectoi@ e^zz and the expressions for the solutions of the second perturbative 



order (Eqs. (4-36), (4-39)-(4-42)) we can rewrite this equation as 



or 



k^q{q - kc 



q{q^ + kl) 



2/>A;2^^M(i)Mi(i) 
1 + /i 



01) 



If we expand the last expression in terms of df'' and keep the lowest order, we find 



^2(1 + 1^) 



m2) 
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The real and the imaginary part have to be satisfied separately and provide the 
scaled growth rate o"^^^ and the scaled frequency uj^^^ as a function of the control 
parameter 



(i)t(i) 



«=0 



(03) 
(|D}4) 

Using the scaled dimensionless growth rate dip = TofJ^^^ with the typical time scale 
To = U2kc{pG + fi2kc)~^, we can rewrite the growth as 

5(i)P(i) ^ _k4iMWM^ 

The choice of the typical time scale r seems arbitrary at this stage, but in ^ where 
we combine the second and the third order, this particular choice is a posteriori 
justified. 

D.2. The contributions in the third order 



05) 



In this section we apply the same arguments as in ^D.l to the normal stress 
boundary condition for the third perturbative order. Taking into account only con- 
tributions proportional to Eq. (B 13) reduces to 

= /Ui/ a^<j(l) _ ^(2)vac^^^(l)vac ^ ^^(1) ^^^(2,1) _ ^(l)vac^^^(2,l)vac 

With the help of the explicit expressions of the eigenfunctions, Eq. ( D[6 ) can be 
written as 



2 ,-. 



2k^ V2{q - kc) ii2 



(2) 



(0)l2 



e«+/,(g^ + A:,^)A2(23^ + 



(2) 



(l)l2^ 



4A;: 



^q^ + kl) 



(2) 



[5, 



(l)l2^ 



(0)l2 



(1) 



^^2+^J'2 



(2) 
t 

(0) 



{l)l2^ 



^2+/i2 dl 



(2) 



[9, 



(0)l2 



(1) 



1 9f) 



V/i2 



P2 + P-2 
1^2 

2_/f_ 
1 + 



[5(0)p 
5(0)l2^ 



e^^^-2A:,>(^2+/i2)::^^ 



«(2) 



5, 



[^(0)13 



(1) 



If expanded in terms of d'f'' we find 

1^2 [a^^'Y 2i/2(l + /i) 



07) 



08) 
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from which uj^"^^ = follows. In the last expression we made use of the results of the 
previous order, namely that u^^^ = in the static limit. For the second contribution 
on the left hand side we can substitute the scaled time derivative d!^^ of the second 
perturbative order and we obtain 

which results in a second order time derivative of the pattern amplitudes. This 
contribution is proportional to the elastic shear modulus fi2 and therefore accounts 
for the reversible bulk processes in the medium whereas the first order time derivative 



represents a purely dissipative process. Equation (D 9) additionally suggests the time 
scale 1^2 1 fJ'2 as the typical time scale to compare oscillatory processes with dissipative 
ones. 

To combine the surface condition with the solvability condition from the bulk 
equations, we multiply by the typical time scale tq and finally obtain (9^^ = tqcj*'^^) 



The second time derivative in Eq. (D 10) is unique for magnetic gels and vanishes in 
the limit of pure ferrofiuids. 
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